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Abstract 

At the LHC, a generic supersymmetric cascade can be a source of top quark. Specifi- 
cally third generation squarks and gluino are the major sources of top quark which could 
also be boosted. In this article, we have shown that jet substructure algorithm can be 
very useful in identifying such boosted top quarks in the cascade. We take inclusive three 
jets plus zero lepton plus missing energy final state and try to reconstruct at least one 
hadronically decaying top quark by using top tagging technique which has good prospect 
at the LHC. 
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I Introduction 



Standard Model (SM) has been successful in explaining a wide variety of phenomena and 
most of the experimental results are consistent with the SM [1] . However, there are some 
issues that make it incomplete as a theory. One of the most important issues is the mass 
of the Higgs boson, the particle responsible for giving mass to the fermions and the gauge 
bosons, is itself radiatively unstable. There is no solution to this puzzle in the SM and 
hence one feels justified to go beyond the SM and look for new physics. Supersymmetry 
(SUSY) [2] is one of the possible extension of SM which can protect Higgs mass and in 
principle solve the problem. Supersymmetry predicts superpartners of SM particles. Since 
no super partner has been observed so far, they must be heavier than the corresponding 
experimental limits j^. This simply implies that the supersymmetry must be broken. 
Unfortunately, the SUSY breaking mechanism is completely unknown. For this reason, 
low energy SUSY spectrum can be completely arbitrary and it needs more than 100 pa- 
rameters to specify the masses and couplings of the SUSY particles. Practically it is very 
difficult to explore such a multi-dimensional parameter space and study its collider signa- 
tures [1]. Therefore, it is necessary to adopt some specific assumptions for SUSY breaking 
mechanisms. There are several phenomenologically viable SUSY breaking mechanisms. 
Gravity mediated supersymmetry breaking model like constrained MSSM (cMSSM) [5] 
is the most popular amongst all. In cMSSM, the full spectrum is completely determined 
only by four and half parameters specified at the high scale. These are universal scalar 
mass (mo), universal gaugino mass (^1/2)1 tri-liner coupling constant (^o); the ratio be- 
tween the vacuum expectation values of up and down type Higgs fields (tan /3) and the 
sign of Higgs mixing term (fi) in the super-potential. This small set of parameters makes 
the model very predictive. 



In a R parit}!^ conserving supersymmetric model, sparticles must be produced in pair 
and eventually cascade down to the lightest SUSY particle (LSP) which can not decay to 
SM particles. Now, if neutralino is the LSP, because of its weakly interacting nature it 
could be a good dark matter candidate. In a collider experiment it is impossible to detect 
such weakly interacting particle, thus contributes as missing energy in the final state. In 
general, the pair production of SUSY particles will lead to final states with multiple jets 
and leptons plus missing energy. So far, many detailed studies have been carried out 
to discover SUSY at the LHC [6], particularly in the context of cMSSM. Such studies 
indicate that inclusive jets plus missing energy channel with zero lepton has the highest 
reach in the cMSSM parameter space [7]. The final state with isolated leptons are also 
important, since the requirement of one or more isolated leptons can suppress huge QCD 
background, although such channels have lower reach than the zero lepton channels [8]. 

If low scale supersymmetry exists in nature, the study of third generation squarks at 
the LHC is of a special interest. This is because, the third generation Yukawa coupling 
is relatively larger than the rest, which results in relatively larger mass splitting between 
the mass eigenstates. Thus the lighter eigenstates could become light enough and might 
be discovered at the Tevatron or in the early stage of LHC [9]. 

parity is a discrete symmetry defined as a Rp~{~lY'^-^^^^^^\ where L, B and S are lepton number, 
baryon number and spin of the particle respectively. 
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Third generation squarks {ii, t2 and 61, 62) can be pair produced at the LHC and 
these can be produced from the decay of SUSY particles. Depending on the parameter 
space, the dominant decay modes of stop squarks (ti, are t Xi 2 3 4 ^^'^ ^ xf2- Simi- 
larly sbottom squarks (61, 62) can decay to t xf 2 ^ Xi 2 3 4- the squarks are heavier 
than gluino except the third generation, gluino will decay entirely to the third generation 
squarks. If all squarks are heavier than gluino, the gluino will decay to three body final 
state through off shell squarks. In general, stop and sbottom squarks, being lighter than 
the first two generation squarks, contribute more to the gluino three body decay. This is 
why, final states with top and bottom quarks in the gluino three body decay are relatively 
more favoured. Final states with b jets accompanying with multiple jets, leptons and 
missing energy have been studied extensively. Another interesting possibility is top rich 
final states from gluino decay. It is thus important to identify the presence of top quark 
in the SUSY cascade, as it carries information about third generation squarks. However, 
top quark identification in a SUSY cascade is not an easy task. Top quark mostly decays 
to W and b quark unless charged Higgs boson is lighter than top quark. In case of lep- 
tonic decay of W, the neutrino, that is present in the final state also contributes to the 
missing energy, in addition to the missing energy from the lightest neutralino. In order to 
reconstruct semileptonically decaying top, we have to separate out the neutrino missing 
energy contribution coming from the top quark. This is why, it is difficult to reconstruct 
a top quark by its semileptonic decay in the cascade. On the other hand, hadronically 
decaying top quark has no real missing energy contribution. However, it is very difficult to 
choose the correct combination of three jets from multiple jets (which are always present 
in a typical SUSY cascade) as it suffers from combinatorial backgrounds. It is thus dif- 
ficult, although, not impossible to identify top quarks in a SUSY cascade. Beside this, 
top quarks in SUSY cascade can be highly boosted due to the usual separation between 
strong and electroweak sector of MSSM. Because of the high boost, the decay products of 
top quark are generally collimated. It makes it very difficult to isolate the decay products 
as separate jets or leptons; and often ends up with a single fat jet in the final state. 

From the above discussion we have seen that the identification of top quark in SUSY 
cascade faces two types of difficulties, e.g., combinatorial background in top reconstruc- 
tion and collimation of the decay products of the highly boosted top quarks. Recently, a 
new method has been proposed to distinguish jets originated from boosted heavy particles 
by the method of jet substructure. By using this method, one can fully reconstruct the 
decaying heavy particle through its hadronic final state. A number of recent studies have 
used this technique to identify boosted top quarks [TOl HH US] , Higgs bosons [l3] , or W/Z 
bosons [U] and the results are very encouraging. It is therefore, a good idea to use the 
above mentioned technique to reconstruct the top signals in the SUSY cascade. 

In this paper we systematically study the top quark production in the SUSY cascade 
and its identification by newly developed top tagging method at the LHC with centre of 
mass energy 14 TeV. We take inclusive three jets plus zero lepton plus missing energy 
final state and try to reconstruct at least one hadronically decaying top quark by using 
top tagging technique, described in [12]. Of course, the production of the boosted top 
will depend on the mass of the parent squark (stop and sbottom), or gluino and their 
corresponding branching fraction to the top quark. In ref [15] authors have studied the 
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pair production of lightest stop squarks assuming the branching ii — )• txi to be equal to 
1 and they have shown that it is possible to extract the stop (ti) - neutralino (x?) mass 
difference. However, this assumption does not hold as we change the parameter space. 

The article is organized as follows. In section II we describe the origin of top quark 
final states from the decay of third generation squarks and gluino and also discuss the 
effective branching to top quark for different scenarios. In section III, we define our bench- 
mark points for the analysis. We discuss our results, briefiy addressing the top tagging 
technique in section IV. Finally we conclude in section V discussing the possible issues of 
this work. 



II Possible Scenarios 

In this section we shall describe the possible sources of top quark from the decay of 
sparticles. In general, the squark mass spectra are very different as one goes from first 
two generations to third generation. The off diagonal term in the squark mass matrix is 
given by mf{Af + fiRf) where Rf = cot /3 for up type squark and Rj = tan/3 for down 
type squark and is the corresponding quark mass |16] . As first two generation quark 
masses are negligible, we can neglect the off diagonal terms in the mass matrix. The 
general hypothesis of flavour blind soft parameter for first two generations avoid poten- 
tially dangerous FCNC and CP-violating effects in MSSM. However, masses of the third 
generation squarks are more complicated because of the Yukawa couplings and the corre- 
sponding changes in RG equations. Due to large mixing through the off diagonal terms, 
there could be large splittings between the mass eigenvalues. Thus, in general the third 
generation masses are quite non-degenerate and lighter ones (ti, 61) can become rather 
light. 

We are now ready to discuss the possible decay modes of third generation squarks. 
If third generation squarks are heavier than gluino, preferably decay to corresponding 
quark and gluino, if kinematically allowed. In case, gluino is heavier than stop or sbottom 
squark, there are several possible interesting channels like 

ii,2 ^ tX? (i=l-4) h,2 ^ bxf (i=l-2) 
h,2 ^ (i=l-4) 61,2 ^ txf (i=l-2) . 

If the mass difference between t2 (62) and ti (bi) or t2 (62) and 61 (ti) is large enough, 
the following decays [T7] 

i2^h ZJHJA h^h W+/H+ 
62 ^ hi Z, /H, I A 62 ^ l\ W+/H+ 

are also possible. Several studies have been performed in this context of CP conserving 
|18] and CP violating MSSM |19j . For a wide range of parameter space, all tree level 
two body decay modes of lightest stop squark are forbidden and the loop decay ti — )■ ex? 
becomes significant [20]. In the case where ti is the NLSP, it decays to three body [21] 
or four body final state [22]. In the next few paragraphs, we discuss the effective decay 
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Figure 1: Branching fraction of ii — )> txi/g for tan (3=10 (left) and tan (3=50 (right) in uiq — mi/2 
plane for Aq = 0, and /i >0. 



branching fractions to top quark from squark and gluino including third generation. We 
use SUSY-HIT |23] for the calculation of decay branching fractions in our analysis. 

ii, 12 to top quark: We have already mentioned the possible decay modes of ii and 12 
squarks. We expect that the top quark contribution from ti will be more than t2, because 
of higher mass value of later. In most of the cases ti is much lighter than t2 and thus ti 
pair production cross section is significant. In order to understand the relevant branching 
of lightest stop squark to top quark, we scan over mQ-mii2 plane of the cMSSM parameter 
space for tan/3 =10 and 50, Aq =0 and /i > 0. We vary mo from 100 to 2000 GeV and 
mij2 from 100 to 1000 GeV. In Figure [H we show the ti to top branching ratio in the mo 
- mi 12 plane. Here we have included ii — )• ty^ and ii — t- tg branching fractions. The white 
region is excluded because fi is the LSP in that region. We have seen that for small mo 
and mi/25 to top quark branching is small (~ 0.2 to 0.3) for both values of tan/3. In 
this region, ti is light and ii — ix? is phase space suppressed and ti — )• bxf is dominant. 
The another reason for it is that, in cMSSM, Xi is often bino like and is wino like 
and bino-squark-quark coupling is small compared to wino-quark-squark coupling. For 
large value of mo and small mi/2) is heavier than gluino and ti — )• tg is dominant due 
to the strong coupling (see yellow region of Figured]). But in this case, stop squarks are 
very heavy and production cross section is negligibly small. For considerably large value 
of stop mass, m^^ is larger than m^o + m^. Since X2 is wino like in most of the cases in 
cMSSM model, ti — ?• X2 ^ is significant (see red region of Figured]). In the similar way, 
t2 squark can also produce top quarks, but, such effects are in general negligible due to 
small cross section. 
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Figure 2: Branching fraction of bi — )■ tx^ for tan/3=10 (left) and tan/3=50 (right) in mo — mi/2 
plane for Aq = 0, and /i >0. 



bi, 62 to top quark: In Figure [21 we show the branching of bi — )■ txi2 the mo - m,i/2 
plane for tan/3 =10 and 50. We can see from the figure that, where ti — )• txi branching 
fraction is large, the branching bi — t- txf is small. This behaviour remains the same 
going from tan/3 = 10 to 50, though for the later case it shrinks a bit. For large niQ 
and small mi/2, higgsino component in the chargino is enhanced because in this region fj, 
becomes small compared to bino (Mi) and wino (M2) soft mass parameters. Otherwise 
as mentioned earlier that Xi is wino like in most of the cMSSM parameter space and for 
this reason bi — t- txf is favoured rather than bi — )• 6x1 ■ This complementary behaviour 
adds to an extra motivation for the search of boosted top quarks in the SUSY cascade in 
which 61 is present. In general, sbottom squarks are heavier than stop squarks and thus, 
the contributions to top quark production may be small. Similarly, t2, &2 also contribute 
to the top quark final state. 

g contribution to top quark: If the gluino mass is above the third generation scalar 
quark mass, gluino decays through g — ?• qq*. The corresponding branching fractions to 
different flavour modes depend on the associated mass values of the quarks and squarks. 
Otherwise if gluino is lighter than squarks, it can still decay to top quark via three-body 
decay [21]. Figure [3] describes the gluino to top branching fraction via two-body and 
three-body decays in cMSSM scenario for tan/3=10 and 50. There are fair amount of 
regions where Br{g — t- ttl) can be very high, whereas, Br(5 — t- bb^) in such regions is very 
low (see Figure [1]) . There is another interesting region in Figure [3| (see conical region) 
where g almost decays to bb^. This is because, in this region, both stop and sbottom 
squarks are lighter than gluino, but, nig — m^^ is less than SM top quark mass. Here the 
only allowed two body decay mode of gluino is bib. The size of such region increases for 
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Figure 3: Branching fraction of ^ to top quark via two body and three body channels for tan/3=10 
(left) and tan (3=50 (right) in the niQ — mi/2 plane for Aq = and /i > 0. 



large value of tan /3 [25] . 

If the gluino mass is lower than the stop and sbottom masses, the gluino decays through 
three body channel to qqx'i or qq xf ■ The ii or hi squarks are generally lighter than other 
squarks and these contribute to the mode ti^ and tbxf. A region of this kind is possible 
for higher values of mo relative to mi/2- In the left black region gluino three body decay 
to top quark is kinematically disallowed. Gluino three body decay to top quarks may play 
a major role in the discovery of focus point region of cMSSM, in which only gauginos are 
light and squarks can be very heavy [26]. 

q contribution to top quark: If squarks are heavier than gluino, these can decay to 
gluino and eventually, gluino decays through three body channels discussed before. In 
this case, the gluino three body final states may contain top quarks. The other possibility 
is that, the gluino can be lighter than first two generations, but, it may be heavier than 
the third generation squarks. In that case, the decay chain like q ^ gq ^ ^i^Q is possible. 

Ill Benchmark Points 

In the previous section, we have discussed the possible sources of top quarks in the SUSY 
cascade in the context of cMSSM. Now, we explore such possibilities by choosing some 
benchmark points in the mg — 1^1/2 plane. The cMSSM model is constrained by different 
low energy observables as well as direct searches like LEP, Tevatron and LHC. We have 
considered such bounds while defining the benchmark points. We scan over niQ -mi/2 
plane for two fixed values of tan/3: 10 and 50 with sign of /U > and ^0 =0. We 
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Figure 4: Branching fraction of g ^ (bbi + c. 
mo — mi/2 plane for Aq = and /i > 0. 
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:) for tan/3=10 (left) and tan/3=50 (right) in the 



take low energy bounds such as 6 — )• 57, Bg — ?• and anomalous magnetic moment 

of muon. We closely follow the reference [27j in our calculation and we use Superlso 
|28] software for this. We also consider recent bounds from 7 TeV LHC data with an 
integrated luminosity 35 pb~^ on cMSSM parameter space in the inclusive multi-jet plus 
missing energy channel. The current experimental lower bound on SM Higgs mass from 
LEP is 114.4 GeV and this bound is applicable to supersymmetric lightest Higgs, except 
for some specific regions in the SUSY parameter space. However, there is an estimated 
theoretical uncertainty of about 3 GeV on the Higgs mass due to higher order effects [29j . 
For this reason we draw separate contours for =111 GeV and 114 GeV. In Figure [5] we 
have shown corresponding bounds discussed above in the mo -mi/2 for tan/3 10 and 50. 
The region bounded by red box is favoured by magnetic moment of muon. The shaded 
region is excluded because in this region fi is the LSP. We have not shown Bg — )• fJ'^fJ'~ 
bound for tan/3=10, because it is much weaker than other bounds considered here. 

We shall tag the top quarks using jet substructure algorithm. This technique can only 
be applicable for highly boosted top quarks. The boost of top quarks in the SUSY cas- 
cade will depend on the relative separation between decaying particle, e.g., gluino, stop 
or sbottom squark and the decay products like neutralino and chargino. The mass differ- 
ence between squark-gluino sector and electroweak gaugino sector is generally large in the 
cMSSM model. To cover different possibilities, we have chosen three specific benchmark 
points for further analysis. 

We choose the following low mass point, which is just above the current LHC bound 



mo = 600 GeV mi/2 = 350GeV Ao = OGeV fi > and tan/3 = 10. 
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Figure 5: Low energy and collider bounds in mo-mi/2 plane for tan/3=10, and tan/3=50. 



For this point the masses of ti, X2 ^^"^ Xi are 675 GeV, 267 GeV and 142 GeV respec- 
tivelyEI- This ensures enough boost to the top, specially the one coming from ii — ?■ txi- In 
this point Br(ti txi)=l8% and Br(ti txl)=W%. This imphes that about 28% of the 
stop squark will decay to a top quark. The stop squark mass is in the higher region and 
the pair production cross section is 93 fb which is much smaller than the total SUSY cross 
section of this point (about 5pb). We also take bibl contribution for the analysis, which is 
about 22 fb, as the parameter point results in slightly heavier sbottom {tti^^ = 846 GeV). 

The branching fraction of bi to the top and chargino mode is ~ 70% (from Figure [2^) . 
The chargino in this case can decay to a and a neutralino and the former can be a 
source of additional jets or leptons. 

Interestingly, there is another mode which could also give rise to top quark. This 
is bi —7- tiW~, where the 6i first decay to ti and the ti decays as usual. As Br(5i — >■ 
iiW~) ~ 8%, contribution to top, neutralino and final state is really small. 

Gluino pair production cross section for this point is ~ 709 fb for nig = 859 GeV, which 
adds to the third generation squarks through its decay, with an extra possible boosted 
top quark. The dominant decay mode of gluino is ^ — )■ tt^ with Br(^ — )■ it^) ~ 95% (see 
Figure [3]). The squark gluino production is about 2.3 pb and it has the major contribution 
to the top final state through the gluino decay to top quark. 

Next, we consider the contribution of gluino decay to sbottom squark. For this purpose, 
we choose the following parameter point from Figure H] as benchmark point 2 

^For our analysis the top is taken to be 173 GeV. 
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mo = 900 GeV mi/2 = 400 GeV Aq = GeV ^ > and tan /3 = 50 



where Br(g — t- bb^ + c.c) = 100%. The hghtest sbottom squark produced from gluino 
further decays to top quark and chargino with branching fraction ~ 59%. This choice 
leads us with rrig = 985 GeV and m^^ = 961 GeV. For this benchmark point, Hghtest stop 
squark mass is about 840 GeV which gives no room for ghiino decays to top quark. 

So far, we have chosen benchmark points which are consistent with low energy con- 
straints and direct search bounds. Now, we consider benchmark point 3 which satisfies 
the dark matter relic density constraint as well as other bounds 

mo = 540 GeV mi/2 = 490 GeV Aq = GeV fj.>0 and tan /3 = 50 . 



The mass of the lightest neutralino is 203 GeV and heavier Higgs masses are about 
480 GeV. The decay widths of heavy Higgs bosons are large and these are about 15-20 
GeV. Neutralinos can annihilate through heavy Higgs mediating s-channel processes and 
satisfy current bound of relic density. We compute the relic density of this point by using 
micrOMEGAs [30] and we find that il/i^ is 0.107 for top quark mass 173 GeV which is 
consistent with WMAP bound nh^ = 0.112 ± 0.007 at 95% CL [31]. 

In this point left handed first two generation squarks are slightly heavier than gluino 
but right handed squarks are lighter. The gluino mass is about 1146 GeV. The branching 
fractions of gluino to stop and sbottom squarks are 36% and 21% respectively. First 
two generation left handed squarks dominantly decay to electroweak gauginos, so, their 
branching to gluino is very small (about 2 — 3%). 

IV Analysis 

The identification of boosted top quarks has been described in the literature. Here we 
briefly describe our algorithm to identify boosted top quarks in the SUSY cascade. We 
have used PYTHIAj33j event generator for generating the events and hadronizations. 
We use FAST JET j34j package for jet formation instead of default subroutine PYCELL 
in PYTHIA. The mass spectrum and decay branching fractions have been generated by 
SUSY-HIT. Here we define set of cuts used in our analysis. 

• We select events with at least 3 jets. The jets are formed by using Cambridge 
Aachen algorithm with a fixed R parameter to be equal to 1 [32] . One may vary this 
parameter to optimize the signal significance. The leading jet px must be greater 
than 300 GeV and other two sub leading jets should have pT greater than 150 and 
100 GeV respectively. Note that, px of the jets are greater than the quoted value 
used by CMS collaboration for SUSY search in inclusive multi-jet plus missing energy 
channels [?]• The absolute value of pseudorapidity of the leading jet must be less 
than 1.7. 
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• The missing pt must be greater than 300 GeV. 

• We put veto on events with isolated leptons with pT greater than 10 GeV. 

• We use variable Ri and R2 which are given by Ri = a/ 6(p2 + (vr — dcpi)"^ and R2 = 
■sj + (tt — 5(f)2Y- Here 5(j)i (502) is the difference between azimuthal angle of 
missing transverse momentum and azimuthal angle of first (second) jet, i.e., J0i,2 = 
l<^/ — 4>j^ 2 1- The i?i and R2 must be greater than 0.5 radian. Also — must 
be greater than 0.3 radian. Also we take \(f^^ — to be greater than 20°. These 
cuts are very useful to reject QCD backgrounds, where the source of missing energy 
is jet energy mis- measurements. 

• The effective mass (M^n) of the event must be greater than 500 GeV. Here Mt is 
defined as the scalar sum of second, third and fourth (if present) jets and 'i/t- 

After passing the above mentioned cuts we choose events for further analysis. We take 
the hardest jet and decluster it into two subjets, if the px of both subjets are greater than 
5% of the parent jet and these are not too close, i.e. [(Jt^I + > 0.1. If the declustering 
is possible, then we take the declustered objects and repeat the same procedure on those 
objects. This procedure is terminated if only one calorimeter cell is left. For our work we 
take a calorimeter cell of 0.1 in both i] and (j). If we start with a top jet (for hadronic top 
decay), we will end up with three or four subjets. The jet mass of the top jet should lie 
in between 150 to 200 GeV. One combination of invariant mass of the two subjets among 
these 3 or 4 jets must be around W boson mass. Also we demand the cosine of helicity 
angle to be greater than 0.7, where the helicity angle is defined as the angle between the 
top quark direction and one of the decay product of W, in the W rest frame. It is also 
possible to tag one subjet as a b jet. However we do not apply b tagging condition in our 
analysis. If the pT of the other sub leading jets are greater than 300 GeV, we apply the 
same algorithm on it and if it satisfies all the criteria, we declare that jet as a top jet. 

In Figure [6] we show the normalized pT distribution of partonic top quarks present in 
the SUSY cascade for three chosen benchmark points. From this plot we have seen that 
a large fraction of top quark events have momenta, greater than 300 GeV. In the right 
hand side of Figure [6] we show the jet mass distribution of events which satisfy all but 
the jet mass cut of the top quark. We observe peaks around top mass in the jet mass 
distributions which verify correctness of the top tagging technique used in the analysis. 
Depending on the parameter space, the peak could be broad as it suffers from SUSY 
contamination in cascade decays. QCD, W/Z+ jets, di-boson + jets backgrounds are not 
taken in the analysis, though the most significant background, ti has been analysed. CMS 
quotes [7] total backgrounds to be of the order of 250 fb for similar type of analysis (with 
less harder cuts used in our analysis) and QCD contribution to it is around 100 fb. Here, 
we have taken much harder cuts on jet pT and missing transverse energy compared to 
CMS cutH and this would further reduce the background. One should note that we have 
used different jet formation algorithm with larger jet cone size. Therefore, the size of the 
SM backgrounds can be different from the CMS analysis. However, the probability of a 
QCD jet, mistagged as a top jet is very small and it is ~ 1% for jet px 300 GeV [12]. 
Further, it is possible to add b-tagging criterion that will also help to reduce QCD and 

^We take leading jet pT to be greater than 300 GeV, where they have used 180 GeV cut on it. For comparison 
see reference [7]. 
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W backgrounds. 

In Table [H we show the individual contributions coming from titl, bibl and gg, qg 
productions in the top signal for the three benchmark points. We do not show the con- 
tributions from t2 and 62 separately, but these are added to the total contribution. We 
may recall that the first two generation squarks are heavier than gluino in our benchmark 
points. These also contribute to the top quark final state through gluino decay. The 
background coming from top quark pair production is not large and it is about 4-5 fb. We 
present our result for integrated luminosity of 30 fb~^. The numbers indicate that even if 
all backgrounds are considered, the signal would be sufficient enough to be observed with 
30 fb~^ of integrated luminosity at the LHC with ^/s = 14 TeV. The Table also shows 
that the gluino contribution is much larger than other contributions like stop and sbottom 
and it is true for all benchmark points. As gluino contribution is too large, it may be 
difficult to determine the hierarchy among the third generation squarks and gluino. A 
detailed study should be worked out to see the prospect of such analysis. 



V Conclusion 

In this paper we have studied the production of boosted top quarks in the SUSY cascade 
at the LHC. In particular, we have taken a specific SUSY breaking mechanism, cMSSM, 
although the outcome of the study is applicable to any general SUSY breaking scenarios or 
any other BSM scenarios. We systematically discussed different possibilities of top quark 
final state in the SUSY cascade. Top quark identification is difficult in cascade decays 
due to the combinatorial backgrounds. Depending on the parameter space, SUSY cascade 
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No. 


Benchmark points 


(30 fb-i) 


(30 fb-i ) 


99 
(30 fb-i) 


9Q 
(30 fb-i) 


Total 
(30 fb-i) 


1 


BPl 


15 


6 


142 


618 


992 


2 


BP2 


8 


3 


110 


336 


591 


3 


BP3 


7 


5 


42 


282 


463 



Table 1: Event rates after top tagging for the benchmark points with an integrated luminosity of 
30 fb~^. The ti contribution is 132 events assuming same integrated luminosity. 



may have top quarks with high p^- We take advantage of this high pT behaviour of top 
quark for its identification by using jet substructure algorithm. This helps to enhance 
the signal significance for the top final state. It has been vindicated in our study that 
for wide range of SUSY parameter space, this technique is fairly applicable. This study 
will help to understand the third generation squarks and gluino masses and decay modes. 
This could also be useful to probe the regions in SUSY parameter space, where squarks 
are inaccessible at the LHC but gluino is in the lighter side. However, we have not done 
detailed background analysis which would be crucial. Also, supersymmetric backgrounds 
may affect the signal significance and it is very much parameter space dependent. One 
should invoke a more sophisticated analysis to overcome such backgrounds which is beyond 
the scope of our study. 

Acknowledgment s : 

The authors would like to thank the organisers of LCWSIO (China, 2010), where initial 
part of this work was done. PB would like to thank KIAS Overseas Travel Grant for 
travel support and BB would like to thank India-ILC forum for travel support to attend 
LCWSIO. BB acknowledge discussions with Sreerup Raychoudhuri, K Sridhar and M. 
Guchait. 

References 

[1] p. Langacker. larXiv:n9ni.0241l [hep-ph]. 

[2] M. Drees, R. Godbole and P. Roy, Hackensack, USA: World Scientific (2004) 555 p; 
For a recent review, see M. E. Peskin, arXiv:0801.1928 [hep-ph]. 

[3] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667, 1 (2008). 

[4] N. Arkani-Hamed, G. L. Kane, J. Thaler and L. T. Wang, JHEP 0608, 070 (2006); 
D. Feldman, Z. Liu and P. Nath, Phys. Rev. Lett. 99, 251802 (2007) [Erratum-ibid. 
100, 069902 (2008)]; C. F. Berger, J. S. Gainer, J. L. Hewett and T. G. Rizzo, JHEP 
0902, 023 (2009); P. Konar, K. T. Matchev, M. Park and G. K. Sarangi, Phys. Rev. 
Lett. 105, 221801 (2010). 

[5] E. Cremmer, B. Julia, J. Scherk, P. van Nieuwenhuizen, S. Ferrara and L. Girardello, 
Phys. Lett. B 79, 231 (1978); E. Cremmer, B. Julia, J. Scherk, S. Ferrara, L. Gi- 
rardello and P. van Nieuwenhuizen, Nucl. Phys. B 147, 105 (1979); R. Barbieri, 



13 



S. Ferrara and C. A. Savoy, Phys. Lett. B 119, 343 (1982); A. H. Chamseddine, 
R. L. Arnowitt and P. Nath, Phys. Rev. Lett. 49, 970 (1982); L. J. Hall, J. D. Lykken 
and S. Weinberg, Phys. Rev. D 27 (1983) 2359; P. Nath, R. L. Arnowitt and 
A. H. Chamseddine, Nucl. Phys. B 227, 121 (1983); N. Ohta, Prog. Theor. Phys. 
70, 542 (1983). 

[6] For recent works, see, for example, S. Bhattacharya, U. Chattopadhyay, D. Choud- 
hury, D. Das and B. Mukhopadhyaya, Phys. Rev. D 81, 075009 (2010); J. Edsjo, 
E. Lundstrom, S. Rydbeck and J. Sjolin, JHEP 1003, 054 (2010); D. Feldman, 
G. Kane, R. Lu and B. D. Nelson, Phys. Lett. B 687, 363 (2010); H. K. Dreiner, 
M. Kramer, J. M. Lindert and B. O'Leary, JHEP 1004, 109 (2010); H. Baer, 
V. Barger, A. Lessa and X. Tata, JHEP 1006, 102 (2010); N. Bhattacharyya, 
A. Datta and S. Poddar, Phys. Rev. D 82, 035003 (2010); U. Chattopadhyay, 
D. Das, D. K. Ghosh and M. Maity, Phys. Rev. D 82, 075013 (2010); ATLAS 
Collaboration, ATL-PHYS-PUB-2010-010 (2010); CMS Collaboration, CMS NOTE 
2010/008 (2010); B. Altunkaynak, M. Holmes, P. Nath, B. D. Nelson and G. Peim, 
larXiv: 1008.3423 [hep-ph] ; B. C. Allanach, S. Grab and H. E. Haber. laDGv: 1010.4261 1 
[hep-ph]. N. Chen, D. Feldman, Z. Liu, P. Nath and G. Peim. larXiv:1011. 12461 [hep- 
ph]. 

[7] T. Yetkin and M. Spiropulu [CMS Collaboration], Acta Phys. Polon. B 38, 661 
(2007); M. Spiropulu, Eur. Phys. J. C 59, 445 (2009); J. B. G. da Costa et al. [Atlas 
Collaboration], Phys. Lett. B 701, 186 (2011) |arXiv:1102.5290] [hep-ex]]; S. Cha- 
trchyan et al. [CMS Collaboration]. larXiv:1107.1279l [hep-ex]. 

[8] Yu. Pakhotin et al, Acta Phys. Polon. B 38, 653 (2007); Yu. Pakhotin et al, CERN- 
CMS-NOTE-2006-134; G. Aad et al [Atlas Collaboration], Phys. Rev. Lett. 106, 
131802 (2011) [ar Xiv: 1102.2357 [hep-ex]]; S. Chatrchyan et al [CMS Collaboration], 
larXiv: 1107. 1 870] [hep-ex]. 

[9] G. Mahlon and G. L. Kane, Phys. Rev. D 55, 2779 (1997); W. Beenakker, M. Kramer, 
T. Plehn, M. Spira and P. M. Zerwas, Nucl. Phys. B 515, 3 (1998); R. Dem- 
ina, J. D. Lykken, K. T. Matchev and A. Nomerotski, Phys. Rev. D 62, 035011 
(2000); A. Dedes and H. K. Dreiner, JHEP 0106, 006 (2001); G. Bozzi, B. Fuks and 
M. Klasen, Phys. Rev. D 72, 035016 (2005); W. Hollik, M. Kollar and M. K. Trenkel, 
JHEP 0802, 018 (2008); M. Beccaria, G. Macorini, L. Panizzi, F. M. Renard and 
C. Verzegnassi, Int. J. Mod. Phys. A 23, 4779 (2008) S. Bornhauser, M. Drees, 
S. Grab et al, [a rXiv :1011.5508 [hep-ph]]. 

[10] K. Agashe, A. B elyaev, T. Krupovnickas et al, Phys. Rev. D77, 015003 (2008). 
| hep-ph/0612015 l. 

[11] M. M. Nojiri and M. Takeuchi, JHEP 0810, 025 (2008); J. Thaler and L. T. Wang, 
JHEP 0807, 092 (2008); L. G. Almeida, S. J. Lee, G. Perez, I. Sung and J. Virzi, 
Phys. Rev. D 79, 074012 (2009); D. Krohn, J. Shelton and L. T. Wang, JHEP 1007, 
041 (2010); B. Bhattacherjee, M. Guchait, S. Raychaudhuri and K. Sridhar, Phys. 
Rev. D 82, 055006 (2010); K. Rehermann and B. Tweedie. la?Xiv:1007.222TI [hep-ph]. 

[12] D. E. Kaplan, K. Rehermann, M. D. Schwartz and B. Tweedie, Phys. Rev. Lett. 101, 
142001 (2008). 

[13] J. M. Butterworth, A. R. Davison, M. Rubin and G. P. Salam, Phys. Rev. Lett. 100, 
242001 (2008); T. Plehn, G. P. Salam, M. Spannowsky, Phys. Rev. Lett. 104, 111801 



14 



(2010); G. D. Kribs, A. Martin, T. S. Roy et al, Phys. Rev. D81, 111501 (2010); 
D. E. Soper, M. Spannowsky, JHEP 1008, 029 (2010); C. -R. Chen, M. M. Nojiri, 
W. Sreethawong, JHEP 1011, 012 (2010); A. Falkowski, D. Krohn, L. T. Wang, 
J. Shelton and A. Thalapillil, larXiv:1006.T650l [hep-ph]; G. D. Kribs, A. Martin, 
T. S. Roy and M. Spannowsky, arXiv: 1006. 1656 [hep-ph]. C. Hackstein and M. Span- 
nowsky, arX iv: 1008.2202 [hep-ph]; J. -H. Kim, [arXiv: 101 1.1493 [hep-ph]]; B. Bel- 
lazzini, C. Csaki, J. Hubisz and J. Shao. [arXiv:1 012.1316 [hep-ph]; 

[14] J. M. Butterworth, J. R. Ehis, A. R. Raklev, JHEP 0705, 033 (2007); L. G. Almeida, 
S. J. Lee, G. Perez et al, Phys. Rev. D79, 074017 (2009); Y. Cui, Z. Han 
and M. D. Schwartz, larXiv: 1012.20771 [hep-ph]; J. Thaler and K. Van Tilburg, 
rarXiv:1011.2268 [hep-ph]. 

[15] T. Plehn, M. Spannowsky, M. Takeuchi and D. Zerwas, JHEP 1010, 078 (2010). 

[16] S. P. Martin, P. Ramond, Phys. Rev. D48 (1993) 5365-5375; S. P. Martin, In *Kane, 
G.L. (ed.): Perspectives on supersymmetry* 1-98. [hep-ph /9709356 1 . 

[17] A. Bartl et al, Phys. Lett. B 435, 118 (1998) ^ arXiv:hep-ph /9804265| . 

[18] A. Datta, M. Guchait, K. K. Jeong, Int. J. Mod. Phys. A14, 2239-2256 (1999); 
A. Djouadi, Y. Mambrini, Phys. Lett. B493, 120-126 (2000); A. Djouadi, M. Guchait, 
Y. Mambrini, Phys. Rev. D64, 095014 (2001); A. Djouadi, Y. Mambrini, Phys. Rev. 
D63, 115005 (2001); A. Belyaev, J. R. Ellis, S. Lola, Phys. Lett. B484, 79-86 (2000); 
A. Datta, A. Djouadi, M. Guchait et al, Phys. Rev. D65, 015007 (2002); M. Graesser, 
J. Shelton, JHEP 0906, 039 (2009). 

[19] P. Langacker, G. Paz, L. -T. Wang et al, JHEP 0707, 055 (2007); A. Bartl, E. Chris- 
tova, K. Hohenwarter-Sodek et al, JHEP 0611, 076 (2006); J. Ellis, F. Moort- 
gat, G. Moortgat-Pick et al, Eur. Phys. J. C60, 633-651 (2009); G. Moortgat-Pick, 
K. Rolbiecki, J. Tattersall et al, JHEP 1001, 004 (2010); G. Moortgat-Pick, K. Rol- 
biecki, J. Tattersall et al, AIP Conf. Proc. 1200, 337-340 (2010); P. Bandyopadhyay, 
|arXiv:1008.3339] [hep-ph]]. 

[20] K. i. Hikasa and M. Kobayashi, Phys. Rev. D 36, 724 (1987). 

[21] W. Porod and T. Wohrmann, Phys. Rev. D 55, 2907 (1997) [Erratum-ibid. D 67, 
059902 (2003)] | :arXiv:hep-ph/9608472| . 

[22] C. Boehm, A. Djouadi, Y. Mambrini, Phys. Rev. D61, 095006 (2000); S. P. Das, 
A. Datta and M. Guchait, Phys. Rev. D 65, 095006 (2002). 

[23] A. Djouadi, M. M. Muhlleitner and M. Spira, Acta Phys. Polon. B 38, 635 (2007) 
[arXiv:h ep-ph/0609 292] . 

[24] H. Baer, V. D. Barger, D. Karatas and X. Tata, Phys. Rev. D 36, 96 (1987); 
R. M. Barnett, J. F. Gunion and H. E. Haber, Phys. Rev. D 37, 1892 (1988); 

A. Bartl, W. Majerotto, B. Mosslacher, N. Oshimo and S. Stippel, Phys. Rev. D 
43, 2214 (1991); J. Hisano, K. Kawagoe, R. Kitano et al, Phys. Rev. D66, 115004 
(2002); P. Gambino, G. F. Giudice, P. Slavich, Nucl. Phys. B726, 35-52 (2005); 

B. S. Acharya, P. Grajek, G. L. Kane, E. Kuflik, K. Suruliz and L. T. Wang, 
V. Barger, W. Y. Keung and B. Yencho, Phys. Lett. B 687, 70 (2010); G. F. Giu- 
dice, T. Han, K. Wang and L. T. Wang, Phys. Rev. D 81, 115011 (2010); A. Bartl, 
K. Hidaka, K. Hohenwarter-Sodek, T. Kernreiter, W. Majerotto and W. Porod, Phys. 



15 



Lett. B 679, 260 (2009) [ariUviOgOS.OTsS' [hep-ph]]; A. Bartl, H. Eberl, E. Ginina, 
B. Herrmann, K. Hidaka, W. Majerotto and W. Porod, arXiv:1107.2775 [hep-ph]. 

[25] A. Bartl, W. Majerotto and W. Porod, Z. Phys. C 64, 499 (1994) [Erratum-ibid. C 
68, 518 (1995)]. 

[26] K. L. Chan, U. Chattopadhyay and P. Nath, Phys. Rev. D 58, 096004 (1998); 
J. L. Feng, K. T. Matchev and T. Moroi, Phys. Rev. D 61, 075005 (2000); U. Chat- 
topadhyay, A. Datta, A. Datta, A. Datta and D. P. Roy, Phys. Lett. B 493, 127 
(2000); J. L. Feng and F. Wilczek, Phys. Lett. B 631, 170 (2005); H. Baer, V. Barger, 
G. Shaughnessy et a/., Phys. Rev. D75, 095010 (2007); S. P. Das, A. Datta, 
M. Guchait, M. Maity and S. Mukherjee, Eur. Phys. J. C 54, 645 (2008); N. Desai 
and B. Mukhopadhyaya, Phys. Rev. D 80, 055019 (2009). 

[27] B. Bhattacherjee, A. Dighe, D. Ghosh and S. Raychaudhuri, Phys. Rev. D 83, 094026 
(2011) [ arXiv:1012.1052l [hep-ph]]. 

[28] F. Mahmoudi, Comput. Phys. Commun. 178, 745 (2008) |arXiv:0710.2T)67] [hep-ph]]: 
F. Mahmoudi, Comput. Phys. Commun. 180, 1579 (2009) [arXiv: 0808 .3144) [heo-ph]] . 

[29] S. Heinemeyer, Int. J. Mod. Phys. A 21, 2659 (2006) |arXiv:hep-ph/0407244 ) . 

[30] G. Belanger, F. Boudjema, A. Pukhov and A. Semenov, Comput. Phys. Commun. 
149, 103 (2002) [arXiv:hep-ph/01122 78|; G. Belanger, F. Boudjema, A. Pukhov and 
A. Semenov, Comput. Phys. Commun. 176, 367 (2007) |arXiv:hep-ph/0607059] . 

[31] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 192, 18 (2011) 
|arXiv:1001.4538] [astro-ph.CO]]. 

[32] Y. L. Dokshitzer, G. D. Leder, S. Moretti and B. R. Webber, JHEP 9708, 001 (1997); 
M. Wobisch, T. Wengler, [h ep-ph/9907280] . 

[33] T. Sjostrand, S. Mrenna and P. Z. Skands, JHEP 0605, 026 (2006) 
[ arXiv:he p-ph /0603175] . 

[34] M. Cacciari and G. P. Salam, Phys. Lett. B 641, 57 (2006) [arXiv:hep-ph/0512210] . 



16 



